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ABSTRACT 

Water masers have been observed in several high redshift active galactic nuclei, in- 
cluding the gravitationally lensed quasar MG 0414+0534. This quasar is lensed into 
four images, and the water maser is detected in two of them. The broadening of the 
maser emission line and its velocity offset are consistent with a group of masers asso- 
ciated with a quasar jet. If the maser group is microlensed we can probe its structure 
and size by observing its microlensing behaviour over time. We present results of a 
high resolution numerical analysis of microlensing of the maser in MG 0414+0534, us- 
ing several physically motivated maser models covering a range of sizes and emission 
profiles. Time- varying spectra of the microlensed maser are generated, displayed, and 
analysed, and the behaviour of the different models compared. The observed maser 
line in MG 0414+0534 is consistent with maser spots as in other quasar jets, provided 
substructure is de-magnified or currently lost in noise; otherwise smooth extended 
maser models are also candidates to generate the observed spectrum. Using measures 
of spectral variability we find that if the maser has small substructure of ~ 0.002 pc 
then a variation of 0.12 mag in flux and 2.0 km s _1 in velocity centroid of the maser 
line could be observed within 2 decades. For the smallest maser model in this study a 
magnification of > 35 is possible 22% of the time, which is of significance in the search 
for other lensed masers. 
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1 INTRODUCTION 

Gravitational lensing occurs when light travelling to 
Earth from distant sources is deflected by the grav- 
itational influence of an int e rvening massive object 
(|Schneider, Ehlers fc Faicol Il992l : IWambsganssI ll99Sl ). The 
light source may be within the Galaxy and lensed by plan- 
etary systems, wh ich is useful in the search for planets (e.g. 
iGaudi et al.ll2008h ; or it may be from a very distant gal axy 



and lensed by anot her galaxy fe.g. IWillis et al.l 120061) or 



galaxy cluster (e.g. ISand. Treu. Elliss et al.ll2008l ). In this 
paper we are concerned with a quasar being lensed by 
a g alaxy, many instances of which have been discovered 
(e.g.lWalsh. Carswell fc Wevmann|[l979l : iHuchra et aTl l 19851; 
iTurner et all 1 19891 : iMvers et all 1 19991 : llnada et all 120051 : 
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iKavo. Inada. Oguri et al1l2010l ). In such systems, multiple 
magnified images of the quasar are produced, with the im- 
age properties determined by the source emission profile, the 
lensing galaxy's mass distribution, and the relative position 
of the source and lens. Assuming a smooth mass distribu- 
tion for the lensing galaxy is sometimes enough to model 
the observed image configurations, but in other cases the ef- 
fect of the compact objects within the galaxy (stars, planets, 
etc.) must also be considered. Due to the intricate nature of 
the light paths through a galaxy of myriad objects, a slight 
change in the location of the background source can produce 
a change in the image magnifications (independent of any 
intrinsic source variability). This effect of the compact struc- 
ture is called microlensing |Chang fc RefsdaJ 19791 : lYound 
ll98ll : |Paczvnskilll986l : rWambsgansslll990l . l2006al ). and identi- 
fiable microlensing fluctuations of high amplitude are called 
high magnification events. These events can las t for rela- 
tively short durations - months or even weeks (|Shalvapinl 



eiy 

I2001D - and usually occur with a frequency of decades. 

The angular size of the source is important for mi- 
crolensing, as larger sources will "wash out" the location- 
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Figure 1. The magnification map representing microlensing variability in image Al, one of the four images of the quasar produced by 
lensing in MG 0414+0534. The map covers a region of the source plane; at all locations the brightness indicates how much a pixel-sized 
source would be (de)magnified if the source were at that location. Bright areas indicate high magnification. There are 6,902,710 Mq 
compact objects in the lens. The map size is 2.5 X 2.5 pc 2 (185 X 185 ER 2 ) with a resolution of 26 AU. The break-out box is 0.1x0.1 
pc 2 (7.4 X 7.4 ER 2 ) at the same resolution. The break-out box shows the high density of caustics in the source plane of this system, 
and that there are patches where there are no caustics and patches where there are many; these appear as light and dark regions on the 
larger map. The microlensing parameters used to generate the map are a (convergence) = 0.472, and 7 (shear) = 0.478 (see Section l2.2ll . 



Table 1. Properties of the magnification maps used. 



Image re 7 Sidclcngth Resolution Number of lens objects (Mq ) 



Al 

A2 



0.472 0.478 



0.485 0.550 



2.5 pc (185 ER) 
15 pc (1115 ER) 
2.5 pc (185 ER) 
15 pc (1115 ER) 



26 AU 
155 AU 
26 AU 
155 AU 



3,204,158 
115,353,274 

6,710,081 
241,562,900 



Four magnification maps were used in this study. Images Al and A2 were modelled, but 
there are two maps of different resolutions for each image because there are a range of sizes 
for the source models, re (convergence) and 7 (shear) are mass parameters from lens models 
for this system, re specifies the effect of mass along a light path, and 7 specifies the effects of 
surrounding mass. 
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induced variability in magnification (|Lewis fc Gil-Merinol 
l2006l ; lBate. Webster fc Wvithj|2008al '), so it is small sources 
that produce the most significant high magnification events. 
The shape of the source emission profile also affect s 
microlensing l|Mortonson. Schechter. fc Wambsganss! 12005?) . 
Therefore, observations and modelling of flux variability can 
provide an estimate of source size and flu x emission profile . 
Apparent chromatic e ffects are also seen (|Lewis et al.lll998l ; 
lEigenbrod et al.ll2008f ) because if an extended source emits 
different frequencies from separate regions, the frequencies 
will be magnified differently, and the source spectrum will 
be altered. Like the magnification, the spectrum may also 
change over time, and studying spectral variability can lead 
to possible spectral emission profiles. In these (and other) 
ways m icrolensing becomes a useful tool for astronomical re- 
search dKavser. Refsdal fc StabelJl986l:IWambsganssll2006bl ; 
iKochanek et al.ll2009l ; IWambsganssll2010? ). 

In this paper we present results of numerical analyses 
of the microlensing of a water maser in a known microlens- 
ing system, MG 0414+0534. We use several physically- 
motivated models to represent the maser source; each model 
is characterised by its shape, size, emission and velocity pro- 
file. A model is then subjected to numerical microlensing so 
that the alteration in the spectrum and flux after lensing can 
be found. The structure of the paper is as follows: Section 
[2] discusses the maser in MG 0414+0534, and mechanisms 
used for numerical microlensing analysis. In Section [3] we 
describe the parameters characterizing the MG 0414+0534 
lens, we specify the maser source models, and explain how 
the models are used. Section [4] presents the results of mi- 
crolensing of the maser by examining which models induce 
the most variability, what time scales may be expected for 
this, what forms the spectra take and how they vary over 
time, and whether the models match the current observa- 
tion. 



2 BACKGROUND 



2.1 The lens and water maser MG 0414+0534 



MG 0414+0534 (|Lawrence. Elston. Januzzi et all Il995l) is 
a gr avitational lens displaying microlensing ( Witt et al.l 
Il995l ). It consists of a quasar at z = 2.639 being lensed 
by a galaxy at z = 0.9584. Recent searches for wa- 
ter masers in six lensed quasars discovered a maser 



in MG 0414+0534 (Impcllizzcri, McK ean. Castangia et al.l 



120081 ; iMcKean et al.l 2010T ). indicating the most distant wa- 
ter found in the universe. The emission comes from the 
616 — 523 transition of H2O at a rest frequency of 22.235 GHz. 
The observations show a line with a single peak, Doppler 
broadened by about ~ 100 km s~\ indicating a spatially 
extended maser source, offset from the systemic velocity of 
the quasar by about -300 km s _1 . There are four images of 
the quasar produced in the MG 0414+0534 system, desig- 
nated Al, A2, B and C, and the maser is observed in Al 
and A2. The sensitivity of the telescope was not enough to 
detect the maser emission in images B and C, which are 
known to be fainter than Al and A2. 

Most known galactic masers are believed to fall into two 
categories: disk or jet masers, both of which may be found 
in the same galaxy and usually consist of groups of "spots" . 



Disk masers are generated in a rotating molecular disk 
around the central engine of an active galaxy and produce 
a characteristic three-peaked spectrum with blue- and red- 
shifted components. T he canonical example of this config- 
urati on is NGC 4258 |Moran. Humphreys. Greenhill et al.l 
2007). Jet masers are probably generated in turbulent shock 
fronts around the edges, or at the ends, of AGN radio jets. 
Such a group of masers will be located some way along 
the jet and be wholly red- or blue-shifted depending on 
the jet direction. The masers are typically offset from the 
AGN core by ~ 100 km s _1 , and have a radial- velocity gra- 
dient of ~ 100 km s _1 mas _1 in a single direction across 
the group. The velocity of the masers is not likely to be 
strongly correlated with the jet velocity for a couple of rea- 
sons. Firstly, the jet speed is likely to be a sig nificant frac- 
tion of the speed of light (|Ros fc K adlcr 200^), and masers 
could not survive in such an environment; secondly, the di- 
rection of the masers may be influenced by turbulent motion 
near t he jet. Examples of galaxies with jet ma sers are NGC 
1052 (IClaussen. Diamond. Braatz et al.Hl998h and Mrk 348 
(|Peck. Henkel. Ulvestad et alj 120031 ). Based on these two 
categories, the maser in MG 0414+0534 is most likely to 
be a jet maser group, since there is only one peak, blue- 
shifted from the core by 300 km s _1 , and a velocity spread 
of ~100 km s" 1 . It could be in a jet that is inclined towards 
us, whi ch is consistent with MG 0414+0534 b eing a type-1 
quasar IjLawrence. Elston. Januzzi et al.lll995h . 

2.2 Analysis of microlensed quasars 

We will adopt a method that has been used previously 
to study microlensi ng in MG 0414+0534 and other 
lense d quasars (e.g. Wambsgans s. Paczvnski fc Schneider] 
19901: ISchechter fc Wambsganss! |2002|: 



Schechter. W ambsg anss fc Lew is 2004; [Lewis fc Gil-Merinol 
20061 : iBate. Floyd. Webster et all l2008bT ). Firstly we cal- 



culate the microlensing of a point-like source (i.e. the 
smallest numerically possible) situated in the source plane 
behind the lensing galaxy. The mass distribution in the 
lens, producing each image, is specified by a convergence 
(k) and a shear (7). The convergence describes the effect 
of mass near a light ray, and the shear incorporates long- 
range influences of the overall mass distribution. Values 
for these parameters for many quasars, incl uding MG 
0414+0534, have been determined pr eviously (|Witt et all 
1 19951 : ISchechter fc Wambsgansi 12002? ) . Based on these, 
point masses are laid down randomly over another grid 
that is the lens plane. Ray-tracing is then used to find the 
path of the light as it travels from the source, is deflected 
through the lens, and reaches the observer. In practice, 
however, it is more computationally efficient to use light 
rays that are fired back from the observer through the 
lens to the source plane. This inverse ray-tracing method 
simultaneously calculates the magnification of the point 
source for many locations over the source plane. 

The result of inverse ray-tracing is a magnification map 
of the source plane, as shown in Figure[T] One map is needed 
for each image. A map does not encode an image of the 
quasar that would be seen by us. Rather, each point on 
the map represents the magnification of a pixel-sized source 
placed at that location, for all locations over the source plane 
region being modelled. Brighter areas are high magnification 
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(a) Grainy Patch Schematic 



T V 

0.002 pc % 



0.13 pc 



0.5 pc 




0.64 pc 



0.53 pc 
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(b) Grainy Ring Schematic 



(c) Hollow & Solid Jets Schematic (1 jet 
shown) 




(d) Grainy Patch: a patch of maser spots (e) Grainy Ring: a ring of maser spots that (f) Hollow Jets: a bi-conical jet structure of 

can be oriented smooth masers that can be oriented. The 

Solid Jet is the same but filled with maser 
material. 

Figure 2. Schematics (a-c) and example images (d-f) of the three complex maser models (the Single Spot and Smooth Patch are not 
shown). Figures (a) and (d) show the Grainy Patch; (b) and (e) the Grainy Ring; (c) and (f) an example of the Hollow Jets, which is 
bi-conical as in (f), but since the jets are mirror images of each other only one schematic is shown in (c). The inclination of the jets in 
(f ) is i = 30 ° South-East, relative to the magnification map as described in the text. The elements of the schematics are not to scale. 



regions and darker areas are low magnification regions, and 
are relative to a mean magnification value. Notice that in the 
map there are thick lines or "trunks" of high magnification 
running up and down, these are the effects of the shear which 
imposes a directional trend on the map. It is not known what 
the true shear direction is for MG 0414+0534. The break- 
out box shows a region of the map at high resolution; one 
can see intricate patterns of light and dark which are called 
caustics. As a source moves across these the magnification 
varies, particularly when the source crosses a bright edge, 
leading to high magnification events. The Einstein Radius 
is defined in the next subsection. 

Spatially extended sources are studied by convolving 
source models with the magnification map . To study the 
differential microlensing dKeeton et al.|[2006l ) of frequencies 
emitted from separate regions of an extended source, each 
region is cut from the source and convolved separately to 
obtain a magnification for that frequency, and then com- 
bined to produce a microlensed spectrum. If this is done for 



points along a path in the source plane then a time-varying 
spectrum for a source is obtained. 

2.3 Microlensing distance and time scales 

Microlensing systems are characterised by several time and 
spatial scales: the time between high magnification events, 
the duration of their occurrence, and the characteristic dis- 
tance scale: the Einstein Radius. An Einstein Ring is gen- 
erated from the simplest lensing situation, where a point 
source is directly behind a point lens, both perfectly in line 
with the observer. In that case the image that is produced 
is a circle around the lens. The observed angle 6e between 
the lens and the circle circumference is given by 

where c the speed of light, G the gravitational constant, 
M the lens mass, Dls,Dl and Ds the distance from lens 
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Figure 3. Measures of variability for the maser line in images Al and A2, for all models, after microlensing. The top panel is the 
variation in the magnitude for image Al, the second panel is the variation in the spectrum velocity centroid for image Al. The next 
two panels are the same measures for image A2. In each panel is a plot of variability measure (vertical axis) as a function of source 
model, and orientation, for those that can be oriented (horizontal axis). Major tick marks indicate the model, inclination direction, and 
inclination angle e.g. GR=Grainy Ring, S=South, ° is the angle. Models that cannot be oriented have no direction or angle. Following 
a major tick mark, minor tick marks indicate increasing inclinations for that model (not all angles shown). The same horizontal axis 
applies to all panels and is only listed on the bottom one. The vertical axis is indicated appropriately for each panel separately. 



to source, observer to lens, and observer to source, respec- 
tively. Projecting 6e onto the source plane produces a dis- 
tance called the Einstein Radius (ER). Thus, the ER incor- 
porates many aspects of the lens system, and is used as a 
length scale in microlensing analysis. It also indicates the 
size scale relevant for sources; sources of order 1 ER in size, 
and smaller, are more lik ely to be susceptibl e to m i crolen sing 
than larger sources (see iRefsdal fc Stabelll (|l99ll . 119971 ) for 
the microlensing of large sources). For MG 0414+0534 we 
calculate the ER based on a point mass of 1 Mq, since that 
size object will be used in our model galaxy. With a redshift 
for the lens of z = 0.9584 and for the quasar of z = 2.639, and 
assuming a concordance cosmology with Qo = 0.3, Ao = 0.7 
and H = 72 km s -1 Mpc -1 , then 1 ER ~ 0.013 pc. The 
ER can be converted to a time scale if the relative veloc- 
ity of the observer, lens, and source is known. Changes over 
time, rather than the ER, represent an observab le of the 
lens system. l|Chartas. Agol. Eracleous et al.ll2002l ) derive a 
velocity of 170 km s _1 for converting source distance to time 
in MG0414+0534, so the time to travel 1 ER is 77 years. 



3 METHOD 



Wc 



the 



nverse ray-shootin g method 
Kavser. Refsdal fc Stabe ll 1986) d eveloped b y Wambsganssl 
1999) and parallelized bv lGarsden fc Lewis] l|2010l ). to pro- 



duce the magnification maps. Source models (detailed 
below) are convolved using Fourier transforms. All of these 
steps are executed using parallel processes and distributed 
memory on a supercomputer. Spectra are then extracted 
from a line across the map chosen to exhibit variability, 
and plotted for graphical examination. Parameters are 
extracted and reduced to measures of variability. 



3.1 Magnification Maps 

The convergence and s hear for MG 0414+0534 are taken 
from IWitt et all l| 19951 ) and listed in Table [U The maser 
has been observed in images Al and A2, so only these will 
be modelled. For each image, two magnification maps are 
used with different resolutions, because the models represent 
masers of various sizes. The highest resolution maps span a 
source plane region of 2.5x2.5 pc 2 (185x 185 ER 2 ), the lower 
resolution maps cover 15 x 15 pc 2 (1115 x 1115 ER 2 ). All 



6 H. Garsden et al. 



maps contain 20000 x 20000 pixels giving a source size reso- 
lution of 26 x 26 AU 2 and 155 x 155 AU 2 respectively, around 
the size of the solar system. All objects comprising the lens 
have masses of 1 solar mass (Mq). The number of objects 
is determined by the lens model and reaches 241,562,900 for 
the low resolution map for image A2. The map shown in Fig- 
ure[T]is the 2.5 x 2.5 pc 2 map for image Al, and the A2 map 
is similar in appearance. We assign directions to the map 
so that the orientation of maser models can be discussed. 
North (N) points up and East (E) is to the right, as in a 
land map. The directions are significant because they relate 
to the direction of shear that produces vertical patterns in 
the maps, which by our convention run N-S. 

3.2 Maser source models 

To model lensing of the maser in MG 0414+0534 we have 
chosen several types of source models of varying levels of 
complexity. We assume the maser material is transparent, 
and emitting isotropically. Because we wish to study the 
effect of microlensing on the spectrum and total flux, each 
model needs a size, a spatial flux emission profile, and a spa- 
tial frequency emission profile. The first two are described in 
Figure The third is specified by the velocity profile of the 
source components, which translates to a frequency emission 
profile due to Doppler shifts. However, the maser velocities 
are not well-understood, as discussed in Section [23] There- 
fore we use the geometry of the model, and estimates based 
on observations, to obtain velocities for the masers, but scale 
them so that the line-of-sight velocity falls in a range of 100 
km s" 1 (except for the Single Spot which has no range). 
This is then sliced into 40 velocity intervals to produce a 
spectrum. All the model descriptions and results will refer 
to the velocities in the model, rather than frequencies. The 
next sections describe the different models. 

3.2.1 Single Spot 

This is the simplest model and is a single maser spot im- 
plemented as a 2-D Gaussian of total width 0.002 pc with 
a central velocity of -300 km s _1 and no velocity broad- 
ening. Since MG 0414+0534 is a powerful quasar we may 
assume that the spot size is at the high end of expected 
sizes for jet masers, so we chose the value of 0.002 pc, 
which also may approximate that of the spots in NGC 1052 
(CI aussen. Diamond. Braatz et al.lll998l ). 

3.2.2 Two patches of masers 

Some masers lie in a group on the side of a jet (e.g. NGC 
1052: IClaussen. Diamond. Braatz et ail Il998l ) with a fairly 
uniform velocity gradient from one side of the group to the 
other. Initially we approach this as an unresolved Smooth 
Patch, modelled as a 2-D Gaussian with a line-of-sight ve- 
locity gradient of -350 to -250 km s" 1 from left to right. 
The total width of the patch is set at 0.1 pc, which is con- 
sistent with the size of other mas er groups (e.g. Mrk348: 
iPeck. Henkel. Ulvestad et all 120031 ). As well as a smooth 
area, we use a group of spots, and replace the Smooth Patch 
with a collection of 49 spots called a Grainy Patch. The ex- 
tent of the patch and the velocity profile remains the same, 



but it consists of Single Spots, randomly located within a cir- 
cular region. Figure [2] (a) has the schematic for the Grainy 
Patch, and (d) shows the resultant model, the Smooth Patch 
is not shown. 



3.2.3 A Ring of Masers 

If the jet in MG 0414+0534 is mostly face-on then it may 
be possible to observe masers situated in a ring around the 
jet. Figure [2] (b) shows the schematic of such a ring, where 
there are 200 Single Spots in a ring of diameter 0.5 pc. The 
location of each spot has been randomly perturbed from 
the true ring. The ring is assumed to lie ~ 1.4 pc along a 
conical jet, with an opening angle of 10 °, allowing the ring 
to be oriented by orienting the jet. The jet has two degrees of 
freedom. Firstly, the angle of inclination determines whether 
the jet is perpendicular or parallel to the plane of the sky, 
in other words parallel to the plane of the magnification 
map. An inclination of i — ° means the jet is pointing 
out of the map towards you, and i = 90 ° is fully side-on. 
Secondly, when the jet is inclined it may be pointing in a 
different compass direction on the map: South/East etc. The 
orientation of the ring has an impact on the lensing, as will 
be shown later. The ring at an inclination of ° is shown in 
Figure rj(e). 

The jet material is outflowing, but the maser material 
in the ring is possibly turbulent and the velocity is not well 
known. However we will assume an overall bulk maser mo- 
tion in a similar direction to the jet. Therefore we set the 
velocity of the masers to be in the same direction as the 
jet material they lie next to, but scaled so the line-of-sight 
velocity of the entire model falls between -350 and -250 km 
s _1 , the same as the for the patches. This range applies re- 
gardless of how the ring is oriented, except when it is face-on, 
in which case the line-of-sight velocity is the same around 
the ring and set to 300 km s _1 . We will call this model the 
Grainy Ring. 



3.2.4 Continuous jet masers 

This model consists of two identical jets in a bi-conical con- 
figuration outflowing from a quasar core. Each jet is 1.59 pc 
long, beginning at a base of width 0.13 pc with an opening 
angle of 25 for 0.53 pc, and then narrowing to 8 for 1.06 
pc, producing a jet diameter of 0.64 pc at the limit of the 
model. The numbers are chosen because of the pixel reso- 
lution of the model. The schematic for one jet is shown in 
Figure [2] (c) . The maser emission region begins 0.66 pc out 
along the jet axis and extends to the end. The maser mate- 
rial can exist all around the jet edge so it forms a truncated 
hollow cone, or it can be distributed through the jet, pro- 
ducing a solid cone. These will be called the "Hollow Jets" 
and "Solid Jets" models. If a cross section is cut though a 
jet perpendicular to its axis producing a ring or disk, the 
flux over that ring or disk remains constant as it travels out 
from the core with the jet material. The bi-conical jet struc- 
ture can be inclined in the same way as the jet carrying the 
ring described in the previous subsection. Since both jets 
are conical outflows we use the same method for obtaining 
the maser velocities as for the Grainy Ring, but with some 
differences. Since the maser material extends along the jets 
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Figure 4. The spectra of the four simplest source models without lcnsing. The y-axis is the flux, the x-axis is the velocity present in the 
source model. Reading left to right, top to bottom, the spectra are for the models: Single Spot, Smooth Patch, Grainy Patch, Grainy 
Ring i = 50° East. The Single Spot has no velocity spread because we have not assigned it any, and produces a sharp peak (within the 
resolution of the slices). The Smooth Patch is a Gaussian and produces a Gaussian spectrum. The patch of spots, each one a Single Spot, 
produces a spectrum of many peaks reflecting the distribution of spots; at the right there is a region where there are a few more spots 
and which produces the highest peak. The ring produces a typical ring-generated spectrum with peaks at the lower and higher ends of 
the spectrum. 



we set the velocity to be decreasing linearly with distance 
outwards, so the fastest material is near the core. Because 
the velocity range is centered around km s _1 with a red- 
and blue-shifted jet on either side, the line-of-sight velocity 
will be scaled to a range of -50 to 50 km s _1 . An example 
of the Hollow Jets inclined to the South-East by i = 30 ° is 
shown in Figure [2] (f ) . 



3.3 Spectrum and Source Path 

We calculate the spectrum as the source travels across a 
region of high magnification on the map, specifically in a 
roughly S-E direction over one of the N-S "trunks" of high 
magnification caustics visible in Figure [T] This means the 
path will begin and end at areas of low magnification but 
pass across caustic regions of high density and complexity 
within the trunk. Such a path is a continuous and representa- 
tive sampling of the microlensing effects that are represented 
by the map. The path is chosen by eye from examining the 
magnification maps at high resolution, and is slightly differ- 
ent for each magnification map used. This method produces 
a time- varying spectrum from which two parameters are ex- 
tracted which are measures of variability: 

(i) A measure of the movement of the velocity centroid of 
the spectrum 

(ii) The variation in total magnitude 

In both cases the standard deviation of the distribution of 
values is used as the variability measure. 



4 RESULTS AND DISCUSSION 
4.1 Time scale 

The magnification map can be used to find the duration 
of, and the time between, high magnification events, which 
occur when the source passes over regions of high magnifi- 
cation in the map. Using a Single Spot for the source, we 
extract a light curve spanning the total width from East to 
West of the map for image Al. Scanning the light curve we 
find instances where the magnitude rises by 1, indicating 
the onset of an event. The shortest distance between two 
such instances is 0.034 pc, giving a time of 195 years be- 
tween events. The shortest distance in which the magnitude 
rises by a value 1 is 0.003 pc, corresponding to a time of 17 
years for the rise time of an event. This translates into sig- 
nificant variability for sources with small substructure like 
spots. For the spectrum of the Grainy Patch model, over a 
period of 17 years, statistical analysis shows that a typical 
variation of 0.12 mag and 2.0 km s _1 in velocity centroid is 
possible. These paths with variability in short intervals are 
more likely to be East- West across the North-South pattern 
of the shear visible in the map. A more detailed analysis is 
beyond this the scope of this paper but could be attempted 
in the future when more informati on is available, and usin g 
more sophisticated techniques (eg. IWvithe fc TurnerlbOOlf ). 



4.2 Model Variability Statistics 

The variability measures indicate how microlensing affects 
a source model as it moves behind the lens, i.e. along the 
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(c) Grainy Patch (d) Grainy Ring, East, 75 

Figure 5. Spectra of the four simplest source models after lensing within image Al, corresponding in the same order to the spectra in 
Figure [4] As the source moves behind the lens a spectrum is captured at each point, these are concatenated to produce a continuous 
spectrum over time, so that the change can be seen. The Velocity axis is the velocity present in the source model. The Flux is the model 
flux. The Time axis indicates the passing of time as the source moves from an arbitrary zero point along a path we have chosen; a Time 
of 600 years corresponds to a source path distance of 0.1 pc (7.8 ER). For each model the spectrum doesn't not change radically; the 
overall shape is retained but small peaks may come and go. The Grainy Patch shows the most variability on statistical measures (Figure 



magnification map on the paths we have chosen. A sum- 
mary of the results is depicted in the panels in Figure [3] 
Each panel is a graph of a variability measure (vertical 
axis) as a function of source model (horizontal axis). The 
top two panels show the light curve and velocity centroid 
variation for all models lensed within image Al, the next 
two are the same, but for image A2. The horizontal axis 
is the same for all panels, and only listed on the bottom 
panel. The vertical axis may be magnitude (mag) or veloc- 
ity (kms -1 ) for each panel, depending on the measure, and 
is labelled on the vertical axis of the panel. The horizon- 
tal axis values are codes indicating the model and orienta- 
tion. Major tick marks are labelled with a code indicating 
the model: SS=Single Spot, SP=Smooth Patch, GR=Grainy 



Ring, SJ=Solid Jet, HJ=Hollow Jet. For the models that can 
be oriented, this code is followed by a direction: S=South, 
S-E=South East, E=East, and then 0°, indicating the ini- 
tial orientation of the model. The minor ticks that follow are 
increasing values of inclination, only 30° and 55 ° are shown. 
The Grainy Ring is inclined by 0°, 10°, 20°, 30°, 30°, 50°, 
55°, 60° and 75°, the Jets by the same, as well as 90°. By 
scanning along the horizontal axis one can immediately see 
the relative variability of all the models, by the height of the 
graph in a panel. 

From these graphs several general comments can be 
made. The Grainy Patch and Ring produce the highest vari- 
ability in velocity centroid, both about the same value of 7 
km s _1 , because they contain spots, i.e. small similar sub- 
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Figure 6. The spectra of the Solid Jets and Hollow Jets models without lensing. Reading left to right, top to bottom, the spectra are 
for the models: Solid Jet i = 0, Solid Jet t = 50 ° South-East, Hollow Jet t = 0, Hollow Jet i = 50 ° South-East. The Solid vs. Hollow 
models produce very similar spectra, because of the way the models are defined. Both of them change when the orientation is changed. A 
face-on (i = 0) inclination (left panels) produces a spectrum with an overall square- wave shape (to the resolution of the velocity slices), 
but when inclined (right panels) this transforms to a single peak, for both models. Only the inclination of i=50 ° South-East is shown 
here as a sample, the peaks get sharper with more inclination. One slight difference is that the Hollow Jet produces a sharper peak than 
the Solid Jets. 



structure. However the Grainy Patch produces more mag- 
nitude variability relative to the Ring because it is smaller 
overall than the Ring. The Single Spot has no velocity cen- 
troid variation because its velocity range is below our res- 
olution, but it has the highest magnitude variability (0.6), 
because it is the smallest-sized model. In contrast, the vari- 
ability for the Jet models is low on both measures, because 
the models are large. The Smooth Patch has variability in 
both measures, less than the Grainy models and more than 
the Jet models, and an unexpectedly high variation in veloc- 
ity centroid, considering it has no substructure. This varia- 
tion will be explained in section [4.3.11 based on its spectral 
behaviour. Using a simple measure of variability per model 
- a sum of the rankings for centroid and magnitude for both 
images - the models in descending order of variability are: 
Grainy Patch, Grainy Ring (best at i=75 ° East), Single 
Spot, Smooth Patch, Hollow Jets, Solid Jets. The overall 
variability is similar for both images Al and A2. 

The orientation relative to the shear changes the vari- 
ability for those models that can be oriented. Any orien- 
tation that produces structures parallel to the direction of 
shear in the map, i.e. N-S, will induce higher microlensing 
variability. This can be seen clearly in the centroid variabil- 
ity produced by the Grainy Ring in image Al (Figure[3] top 
panel). As the direction changes from South to East, and 
with a high inclination, the ring presents a thin N-S pro- 
file, so the microlensing variability increases. In the case of 
the Jets, when they are N-S and highly inclined, the entire 
structure is in the direction of shear and so shows the most 
magnitude variability. 



4.3 The forms of the spectra 

The jets will be discussed separately from the other models 
since they have a different velocity spread and we will see 
they behave differently under inclinations. 

4-3.1 Spot, Patches and Ring 

The form of the spectra for these source models is given in 
Figure [4] The Single Spot shows a single peak with no ve- 
locity spread. The Smooth Patch produces a Gaussian spec- 
trum, as it is a Gaussian model. The Grainy Patch has a 
multi-peaked spectrum due to the many small spots. The 
Grainy Ring shows the typical spectrum shape generated by 
a ring structure, with two peaks on the sides of the spectrum 
and a continuum between them, it is slightly broken because 
of the grainy nature. The spectra in Figure [4] are good repre- 
sentative examples as the spectra for all models retain those 
overall shapes even when inclined (remember the velocities 
have been normalised to always lie in the range of -350 to 
-250 km s _1 ), and they also retain that overall shape when 
lensed. The changes that occur due to microlensing can be 
described either as an overall rise and fall or the introduc- 
tion/disappearance/movement/magnification of small peaks 
or folds. Figure [5] contain typical examples of these phenom- 
ena, each spectrum is the lensed spectrum over time corre- 
sponding to the unlensed spectrum of Figure [4] The Single 
Spot shows how the magnification can rise and drop signif- 
icantly when the quasar is travelling across a trunk of high 
magnification on the magnification map. The peak in the 
velocity axis rises from Time = yrs to a maximum when 
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(c) Hollow Jets, i = (d) Hollow Jets, i = 50 ° South-East 



Figure 7. Spectra of the Jets source models, for image Al, after lensing. Each image corresponds to the spectrum in Figure|5] As the 
source moves behind the lens a spectrum is captured at each point, these arc concatenated to produce a continuous spectrum over time, 
allowing the changes to be seen. The Velocity axis is the velocity present in the source model. The Flux is the model flux. The Time 
axis indicates the passing of time as the source moves from an arbitrary point along a path we have chosen; a Time of 1500 years 
corresponds to a source path distance of 0.26 pc (19 ER). The models produce different spectra at different orientations but these are 
little affected by microlensing, and change little over time. 



Time ~ 500 yrs, and then drops sharply after that. The 
Smooth Patch when lensed has a similar peak that rises and 
falls over time, but in the Velocity axis the peak is wide and 
smooth. One can also see that at Time=0 yrs, Velocity ~ - 
280 km s~ x there is a fold, which diminishes over time. Also, 
although it is difficult to see in Figure [SJb), the entire peak 
for the Smooth Patch drifts from left to right on the Velocity 
axis, explaining the high centroid variability for this model; 
otherwise, it would have a very low centroid variability. The 
Grainy Patch has many peaks which rise and fall at differ- 
ent times, due to the substructure in the model, with an 
overall rise and fall as in the previous spectra. The Grainy 
Ring varies in flux but only shows some rise and fall at the 
edges, and at different times; it also shows an interior ridge 
formed by a peak that appears near Velocity ~ -260 km s _1 



at Time ~ 400 yrs and then moves to the other side of the 
spectrum arriving at Velocity ~ -340 km s _1 at Time ~ 700 
yrs, mimicking the drift in the Smooth Patch spectra. Other 
Ring orientations manifest similar behaviours. The rise and 
fall of the spectra, as in the Single Spot, will affect the total 
observed brightness of the maser line over time; the change 
of folds and peaks will affect the location of the velocity cen- 
troid of the broadened maser line over time. We do not see 
any dramatic changes in spectrum shape, based on viewing 
all the spectra obtained. 

4.3.2 Jets 

The jets are a bi-conical structure with identical outflowing 
jets on either side. The velocity profile is a linear veloc- 
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ity gradient out from the core, red- or blue-shifted on each 
side. With this model the form of the unlensed spectrum 
can change significantly depending on the orientation, but 
not over time, in contrast to the other models. Figure [6] 
shows two Solid Jets models and two Hollow Jets models 
before lensing. The Solid Jets and Hollow Jets when face-on 
produce an identical spectrum because the fluxes are nor- 
malized to be the same, and in this orientation the observed 
velocities will be the same. The apparent peaks are due to 
the resolution of the velocity slicing, if the resolution was 
very small the spectra would be smooth and flat. When the 
same models are inclined by 50 ° , they both transform into a 
single peak, which becomes sharper the more inclined they 
are. Figure [7] shows the same models after lensing, and it 
is clear there is little change due to lensing, because the jet 
structure is large compared to the Einstein Radius and does 
not undergo significant microlensing. There is some varia- 
tion right at the peak with small folds coming or going, 
but the flux is fairly constant, indicating little variation in 
the spectrum over time. What is of more interest with this 
structure is the spectra produced by different inclinations. 
Starting with a square wave when face-on, both the Solid 
and Hollow Jets spectra can merge towards a single smooth 
peak whether lensed or unlensed, with microlensing produc- 
ing slight variations at the peak over time. This has implica- 
tions for matching source models with the observed spectra. 

4.4 Comparison with observations 

Several factors make it difficult to choose which of our mod- 
els can produce the observed maser line. Firstly, in a source 
with multiple components, some components may be de- 
magnified by microlensing and thus not appear in the spec- 
trum; observation over long time scales may reveal these. 
Secondly, the entire maser spectrum may not be clearly vis- 
ible because of noise in the signal or components of the 
maser spectrum may lie below the detection limit. Thirdly, 
the observed line is coincident with an HI absorption trough 
IIBarvainis et al II 19981 ; I ICurran et alj2007h , possibly masking 
a different spectral signature. Fourthly, the observed line is 
the combined spectra of both image Al and A2. Although 
the images have similar microlensing properties and simi- 
lar variability measures, the spatial correlation between the 
high resolution maps for Al and A2 is only 56%. This indi- 
cates that some high magnification events do not correlate 
between the two images in either magnification or epoch of 
observation. It would be preferable to observe the maser line 
in Al and A2 separately, since this will provide extra infor- 
mation to use to distinguish between models. Naturally the 
last three issues are due to the technical limitations of the 
observation, and will improve in the future. 

To match the observation, we look for models that have 
produced a single peak. The Single Spot has a peak, but it is 
not likely to be realistic since masers in jets usually occur in 
groups of spots. The Smooth Patch produces a smooth peak 
that is broad and maintains its shape, rising and falling with 
changing magnification, although it may drift around in ve- 
locity. Therefore if observation shows little change in shape, 
but changes in flux and centroid, the Smooth Patch is indi- 
cated. The Grainy Patch produces several peaks depending 
on the arrangement of maser spots, and our model deliber- 
ately includes a large peak which could match the observa- 



tion. The existence of other smaller peaks in the model spec- 
trum, but not seen in the observed spectrum, does not lessen 
the likelihood for the Grainy Patch. Such peaks may be cur- 
rently de-magnified or lie below the noise, either within the 
current observed line or at lower or higher velocities outside 
it. In the latter case the Grainy Patch model remains valid 
as our velocity range is artificial and can also be expanded. If 
future observations resolve smaller peaks, or peaks appear 
over time, the Grainy Patch model is consistent with the 
observation. The Grainy Ring produces two clear peaks at 
the higher and lower ends of the spectrum and we do not 
observe enough microlensing variation that could produce a 
single peak, so this model appears less likely. Interestingly, 
the Jet models are capable of producing a smooth single 
peak, which varies insignificantly over time in either velocity 
or flux. If the observed maser spectrum exhibited the same 
properties, then the Jet models would be consistent with the 
observation. However, an inclination of more than 50 ° is re- 
quired to produce such a spectrum, and since the quasar is 
believed to be not so inclined, this mitigates against the Jet 
models. 



4.5 Magnifications 

A search f or masers in other lensed AGNs has not been 
successful (|McKean et al.l [20101 ). partly due to the lack of 
required sensitivity in existing telescopes, but that will im- 
prove. However, predictions have been made about the de- 
tectability of lensed AGN masers assuming a magnification 
of 35 for M G 0414+0534 based on analytic strong lens- 
ing models |lmpellizzeri. McKean. Castangia et al.l 120081 : 
iTrotter. Winn fc Hewittl 120001 ) .Using the Single Spot as a 
source, we found a peak magnification of 91 can be achieved 
in image A2, with a magnification of 35 or better achieved 
22% of the time. If long-term observing, combined with mod- 
els of the system, show that the typical magnification is 
lower than 35, this indicates the maser is more intrinsically 
luminous than previously thought, with consequences for the 
rate of detection. 



5 CONCLUSIONS 

Observations of jet masers in active galaxies indicate that 
the masers typically exist in a group with a velocity range 
of the order 100 km s _1 . The maser in MG 0414+0534 con- 
sists of a single peak. This is not consistent with a group of 
maser spots that produces multiple peaks, unless some peaks 
are lost in noise, or currently de-magnified, which is cer- 
tainly possible. On the other hand the observation could be 
a smooth patch of masers, or a smooth distribution through 
the inner parts of the quasar jets, all of which also produce 
a single peak. A patch of spots will produce the most vari- 
ability, in both velocity centroid and magnitude, the smooth 
patch will produce less - it may drift in velocity but remain 
fairly constant in shape - the jets will not change at all. A 
source with substructure like a patch of spots could produce 
changes in the spectrum within a couple of hundred years, 
lasting for 10-20 years. These results are just the beginning 
of the numerical analysis of this lensed maser, and continu- 
ing work would include: 
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• Developing physically consistent models based on fu- 
ture observations and new understanding of masers 

• Developing statistical tools for the extraction of his- 
tograms, structure functions, and other data from the high 
resolution magnification maps 

but it would be prudent to wait for observations that have a 
higher sensitivity and resolution of flux and velocity, and in 
particular that are able to detect the maser in images B and 
C. A long-term monitoring program should be carried out 
on this source, either with a large single dish or a sensitive 
VLBI array (e.g. the European VLBI Network) with higher 
resolution than the original Effelsberg observations, with the 
following aims: 

• monitoring should be conducted on monthly time scales 
to pick up short term fluctuations as well as building statis- 
tics for long term variability 

• the Al and A2 image spectra should be separated, 
which allow a comparison of the spectral behaviour between 
the two images for all the models, adding information that 
may constrain the choice of model 

• observations at higher spectral resolution may allow a 
more accurate determination of the spectral shape; differen- 
tial fluctuations over time due to substructure will help to 
determine the morphology of a maser group. 
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